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We investigate the dynamics of quantum and classical correlations in a system of two 
qubits under local colored- noise dephasing channels. The time evolution of a single qubit 
interacting v^ith its own environment is described by a memory kernel non-Markovian 
master equation. The memory effects of the non-Markovian reservoirs introduce new 
features in the dynamics of quantum and classical correlations compared to the white 
noise Markovian case. Depending on the geometry of the initial state the system can 
exhibit frozen discord and multiple sudden transitions between classical and quantum 
decoherence [L. Mazzola, J. Piilo, S. Maniscalco, Phys. Rev. Lett. 104, 200401 (2010)]. 
We provide a geometric interpretation of those phenomena in terms of the distance of 
the state under investigation to its closest classical state in the Hilbert space of the 
system. 

Keywords: discord; classical correlation; non-Markovian. 

1. Introduction 

In the last two decades lots of interest have been devoted to the definition and 
understanding of correlations in quantum systems. From the first formalization 
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of the separability problem by Werner ^ a number of fundamental results have 
been found in this field. In particular recently a big deal of attention has been 
devoted to the definition and study of quantum and classical correlations in quan- 
tum systems. The works of Ollivier and Zurek and Henderson and Vedral ^ 
had marked the beginning of a new line of research shifting the attention from the 
entanglement vs. separability dichotomy to the quantum vs. classical paradigm. 
The fact that separable states can be prepared via local operations and classical 
communications between two parties has often been interpreted as a signature of 
classicality. However, there are quantum correlations which are not captured by en- 
tanglement. For example the state p = 1/2(|0)(0|^ (g) + |+)(+|^«) |1)(1|b) 
with |±) = (|0) ± \1))/V2 is separable, nevertheless it cannot be described by clas- 
sical means, i.e., by a classical bivariative probability distribution. The reason for 
that is hidden in the non-orthogonality of the states of A and B, and consequently 
in the impossibility to locally distinguish the states of each subsystem. 

Many different measures have been proposed to describe such more-general- 
than-entanglement quantum correlations-^ - , among them probably the most 
popular one is the quantum discord The discord, defined as the difference be- 
tween the quantum generalization of two equivalent formulations of classical mutual 
information, involves an optimization procedure over the set of measurements on a 
given subsystem. Therefore mathematical investigations have been carried to find 
analytical expressions fo r the d iscord in different types of systems, as in qubits^^^ 
and harmonic oscillators ^^^^l^^. The set of zero discord states has been shown to be 
measure zero and nowhere dense moreover a new class of multipartite separable 
states, the pseud o-cl assical states, has been introduced and showed to have inter- 
esting propertie s P roposals to witness the presence of quantum discord have 
appeared in Ref. I ^'^ ^ ^^ i The role of quantum discord has been investigated in quan- 
tum information tasks as the power-of-one-qubit protocol!^ or in the Grover search 
algorithm posing the possibility of using quantum discord as a new, more gen- 
eral than entanglement, quantum computation resource. The dynamics of quantum 
discord has been studied even in quantum biology, specifically in light-harvesting 
complex 

A very active area of investigation of the discord is its behaviour under de- 
coherence, i.e . the dynamics of qu antum and classical correlations in op en quan- 
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tum systems I 19 | 20 | 21 | 22 | 23 | 24 | 25r2 6 27 28 29 30 particular in Ref. ED we have 
discovered that quantum correlations, quantified by quantum discord, can be com- 
pletely unaffected by decoherence for long intervals of time. Specifically, we have 
found that an open system of two qubits, in Markovian dephasing channels can 
undergo, what we have called a sudden transition between classical and quantum 
decoherence, namely a transition between a "classical decoherence" phase in which 
quantum correlation is frozen, while classical correlation is lost, and a "quantum 
decoherence" regime in which quantum correlation is deteriorated while classical 
correlation remains constant. 
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Here we study the time evolution of classical and quantum correlations for a 
system of two identical qubits under local non-dissipative non-Markovian channels. 
We use a well established model describing the dynamics of a two-leve l sy stem in 
a dephasing channel in which white noise is replaced by colored noise There, 
the interaction with the environment is schematized by random telegraph signal 
noise. We study the non-trivial effects the non-Markovianity of the system brings 
to light. As in the Markovian case, there exists a class of initial states for which the 
system exhibits the su dden tr ansition between classical and quantum decoherence 
and the frozen discord I 3f I 30 l However, here multiple transitions can arise due to 
the memory effects of the environment. 

The transition between the two dynamical regimes can be explained, in the light 
of the results of Modi et al. There, the discord is defined as the relative entropy 
with respect to the closest classical state, therefore frozen discord corresponds to 
a dynamical evolution in which the state under investigation maintains a fixed 
distance (relative entropy or trace distance) to its closest classical state. In this 
framework the sudden transition corresponds to a sudden change in the location of 
the closest classical state in the Hilbert space of the system. 

The outline of the paper is the following. In Sec. II we introduce the physical 
system under investigation I^H. In Sec. HI we recall the concepts of quantum dis- 
cord and classical correlation and study the dynamics of those quantities in our 
specific system. In Sec. IV we illustrate the interpretation of the phenomenon of 
frozen discord in terms of the distance to the closest classical state set. Section V 
summarizes and concludes the paper. 

2. The model 

We begin by introducing the physical model under study: two non-interacting qubits 
under local identical colored noise dephasing channels. Since the two qubits are 
locally interacting with their own environments, and are not coupled between each 
other, they have independent dynamical evolutions. Therefore we can first consider 
the dynamics of the single qubit, and then derive the dynamics of the composite 
two-qubit system. 

A generic memory kernel master equation can be written in the form 

P^KCp, (1) 

where p is the density operator of the small system of interest, in our case a qubit, 
£ is a Lindblad superoperator describing the dynamics induced by the interaction 
with the environment, K is an integral operator acting in the following way K<j) — 
Jq k{t — t')(f>(t')dt' , and k{t — t') is the kernel function determining the type of 
memory of the reservoir. 

A master equation of this form arises when considering a two-level system sub- 
jected to random telegraphic noise. A model like this describes, for example, a 
spin in presence of a magnetic field, having constant intensity but inverting its 
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sign randomly in time. It is possible to write a time-dependent phenomenological 
Hamiltonian for this kind of system: 

3 

H{t) = hY^T,it)a, (2) 

i=l 

here ai are the Pauli operators and Ti{t) are independent random variables. We 
consider the general case of noise in the three directions. Each random variable can 
be written as Ti{t) = aiJiiit). The random variable ni{t) has a Poisson distribution 
with a mean equal to t/2Ti, while is a coin- flip random variable assuming the 
values ±ai. By considering the von Neumann equation p = —{i/h)[H, p], one can 
find a formal solution for the qubit density matrix operator of the form 

Pit) = p(0) - z / Y.^kis)[ak,pis)]ds. (3) 

When inserting such a formal solution back into the von Neumann equation and 
performing a stochastic average, one obtains the following memory kernel master 
equation 

P(t) = - f ^exp(-(i-t')/Tfe)a2[afe,[afe,p(t')]rfi' (4) 
•^0 k 

in which the memory kernel comes from the correlation functions of the random 
telegraph signals (F^ (t)rfc(t')) = al cxp {-\t - t'\/Tk) S^k- 

The dynamics of this system and in particular the conditions of complete pos- 
itivity of the map corresponding to such a master equation have been studied in 
detail by Daffer et al. in Ref. In particular they demonstrated that complete 
positivity is assured when two of the are zero, representing the physical situation 
of noise in only one direction. This is the case we are focusing here. In that work 
they also provided a Kraus operator representation of the map ^t{p) = J2k ^kP^k- 
Depending on which of the = a is non zero the map has Kraus operators 

A, =0, Ak = 0, (5) 
Ai =^[1 + A{iy)]/2I, 

where A(i') — exp zy[cos (/ii^) + sinipiy)/ p], p = ^ (Aar)^ — 1, u ~ i/2r is a dimen- 
sionless time, i = 1, 2, 3 is the direction of noise, and j and k the directions in 
which there is no noise. So by changing the direction of the noise we get a colored 
noise bit flip, bit-phase flip or phase flip channel, respectively. Having the Kraus 
operators of a single qubit, one can write the evolution of the generic state of the 
two qubits each interacting locally with its own environment as 
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where is the completely positive trace preserving map ruling the evolution of 
the system. 

It is worth mentioning that, even though in general the mathematical structure 
of the memory kernel master equations m E g. (j4|) do not guarantee the presence of 
non-Markovian features in the dynamics'^, however the quantum process descri bed 
by $t(-) is non-Markovian according to the measure of non-Markovianity of Ref.^^. 



3. Dynamics of quantum and classical correlations 

The physical quantities we want to investigate here are the quantum and classical 
correlations between the two qubits. We use the quantum discord as a measure of 
quantum correlations. 

The idea behind quantum discord is to exploit the difference between quan- 
tum extensions of classically equivalent concepts to evaluate the "quantumness" of 
a quantum system. In particular in classical information theory there are two equiv- 
alent expressions for the mutual information of a bipartite system. However, when 
pursuing the quantum analogue, these two formulations differ, and one can use the 
mismatch between the two quantum extensions of classical mutual information to 
assess quantum correlations. 

The first quantum generalization of mutual information is the so called quantum 
mutual information: 

I{pab) = S{pa) + S{pb) - S{pab), (7) 

where pab is the density matrix of the total system, Pa(b) is the reduced density 
matrix of subsystem A{B)^ and S{p) = — Tr{plog2 p} is the von Neumann entropy. 
This is generally accepted as the measure for the total amount of correlations (quan- 
tum and classical) of a quantum system ^ . 

The second extension of mutual information requires the generalization of the 
concept of conditional entropy. Indeed, performing measurements on system B af- 
fects our knowledge of system A, in particular how much system A is modified by 
the measurement depends on the choice of the measurement performed on B. Here 
the measurement Bk is considered of von Neumann type and it is described by a 
complete set of orthonormal projectors H^ corresponding to outcome k. So the con- 
ditional density operator, which is the quantum state of the total system conditioned 
on the measurement outcome labeled by k, becomes Pk — {I ® Bk)pAB [I ® Bk)/pk 
where pfc — Tr{(I(X)Bk)/CAB(I'8)Bk)}, and / is the identity operator for subsystem A. 
Defining the quantum analog of conditional entropy as S{pAB\{Bk}) = '^kPk^(Pk) 
we can introduce the second quantum extension of mutual information. 

JipAB\Bk) = S{pa) - SipAB\Bk). (8) 

In Ref . El Henderson and Vedral have shown that the maximum of this quantity over 
all the possible set of measurements can be interpreted as a measure of classical 
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correlations of the state 

C{pab) = ma.x[J{pAB\Bk)]. (9) 

Therefore, the difference between the the quantum mutual information I{pab), 
describing total correlations, and the measurement-based definition of quantum 
mutual information C{pab), measuring classical correlations, defines the so-called 
quantum discord 

D{pAB)^IipAB)~C{pAB)- (10) 

In the following we investigate the dynamics of these quantities in our non- 
Markovian dynamical model for a particular class of states. The states we consider 
have maximally mixed marginals and have the form 



PAB = 7 ( Iab + 1 ' (11) 



where crf^^'^ are the Pauli operators in direction i acting on A (B), Ci are real 
numbers such that < |ci| < 1 for every i, and Iab the identity operator of system 
AB. This class of states can equivalently be written in the form of Bell diagonal 
states 

PAB{t) = A+(t)|vl/ + )(vl,+ | + A+(t)|$+)($+| 

+ X-{t)\^-){<i>-\+\-{t)\^-){^~\, (12) 

where 

A± (t) = [1 ± ci it) T C2 [t) + C3 m /4, (13) 
A± (t) - [1 ± ci (t) ± C2 [t) - C3 (i )] /4, (14) 

and 1^-=^) = (|00) ± |ll))/y2, = (|01) ± |10))/y2 arc the four Bell states. 

The expressions of the classical and quantum correlations for this class of states 
are given by^^ 

C{pab) = E ^-^^Y^ ^°S2(1 + (-l)'x) (15) 

/c=l 

D{pab) = 2 + ^ \\{t) log2 \\{t) - C{pab), (16) 

k,l 

where x{i) = max{|ci(t)|, \c2{t)\, |c3(t)|}, fc = $, and I = ±. 

Since the dynamical evolution map $t(') of Eq. © docs not change the form 
of the state in Eq. (|lip . classical and quantum correlations are given by those 
expressions throughout all the time evolution. 

As already noticed at the end of Sec. II, the interaction with the environment 
gives rise to a bit flip, bit-phase flip or phase flip channel when the direction of 
the noise is along x, y or z (in Eq. JS]) i — 1,2,3), respectively. Using the Kraus 
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operators in Eq. ([5]) and the dynamical evolution of the state in Eq. (j6]) we find 
that the parameters Ci{t) evolve in the following way 

c^{u) = c,(0), Cj{u) = Cj{Q)K{uf, Ck{v) = Cfe(0)A(iy)2, (17) 

where v — ijlT is dimensionless time, i is the direction of the noise, and j and fc the 
other two directions. We focus on the phase flip channel (which means z = 3) case 
since completely equivalent results are found for the other two types of channels. 
The corresponding results for bit flip (bit-phase flip) channel are simply found 
exchanging C3 with c\ (c-i). 

As in the Markovian case one can distinguish three difi^erent dynamical 

regimes, depending on the relations between ci, C2, and C3. 

First regime: |c3(0)| > |ci(0)|, |c2(0)|. In this regime, classical correlations re- 
main constant in time. This is clearly displayed by Eq. (llSp . where x(i) is equal to 
1 03(0) I which is constant during all the time evolution. The quantum mutual infor- 
mation characterizing the total correlations displays damped oscillations, therefore 
quantum discord is also exhibiting damped oscillations, tending asymptotically to 
zero. No discontinuities appear in quantum or classical correlations dynamics. 

Second regime: C3 = 0. Here, all the three types of correlations (quantum, clas- 
sical, total) display damped oscillations and tend asymptotically to zero, meaning 
that the state of the system tends asymptotically to a product state. 

Third regime: |c3(0)| < |ci(0)| and/or |c2(0)|. This is the region displaying most 
interesting dynamical features. Initially, classical corre lati ons decay till to a cer- 
tain point in time, and then become abruptly constant and at that very same 
time the quantum discord changes its decay rate. In Ref. we demonstrated that 
the initial decay rate can be even equal to zero. This is mathematically explained 
observing the analytic formula of classical correlation in Eq. ()15|) : initially x(i) is 
equal to max{|ci(i)|, |c2(i)|} with ci{t) and C2{t) oscillating functions, then at the 
sudden change point |c3(0)| becomes bigger than max{|ci(i)|, |c2(t)|}, therefore x{t) 
is now equal to |c3(0)| which is constant in time, hence classical correlation becomes 
constant. 

As time passes new features may appear in the dynamics of correlations. In 
fact if during the time evolution |ci(t)| or |c2(t)| turns bigger than |c3(0)| classical 
correlation is not constant anymore, and starts again oscillating until it becomes 
again abruptly constant. Clearly, sudden changes in the dynamics of discord take 
place in correspondence of the discontinuities of x(t). Within this regime we can 
distinguish two different behaviours for the dynamics of the quantum correlation: 
the frozen behaviour, and the sudden change dynamics. 

(i) Frozen discord with sudden changes displayed in Fig. [T] (a). The dynamics 
displays multiple sudden transitions: classical correlations at flrst decrease, while 
the discord remains constant unaffected by decoherence, then at a sudden transition 
point classical correlation becomes constant and discord oscillates, and then again 
discord becomes constant while classical correlation oscillates. The initial state in 
Fig.[T](a) is of the form pab = (1 + C3)/2|*+)(*+| + (1 - C3)/2|$+)($+| with C3 = 
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Fig. 1. (Color online) Dynamics of mutual information (green dotted line), classical correlations 
(red dashed line) and quantum discord (blue solid line) of two qubits locally interacting with non- 
Mar kovian dephasing channels as a function of scaled time u = t/2T with r = 5 s and a = 1 . 
(a) Frozen discord and multiple transition regime, here we have set ci(0) = 1, C2(0) = —0.6 and 
C3 = 0.6. (b) Sudden change regime, obtained setting ci(0) = 0.35, C2(0) = —0.3 and C3 = 0.1; 
the arrows emphasize the sudden change points. 
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0.6, however a more general class of states exhibits frozen discord and (multiple) 
sudden transitions. This class of states is of the form of Eq. (fTTj) where ci(2)(0) = k, 
C2(i)(0) = — cafc and 03(0) = C3 with k real and |fc| > |c3|. For this class of states 
the quantum mutual information describing the total correlation of the system can 
be cast in a very simple form: 

2 

ApABm = ^'^^s^^'"' l0g2[l + (18) 

^ l + (-l)^Ci(2)(t) , 
J = l 

where the first constant term represents either the frozen discord or the constant 
classical correlation depending on the dynamical phase of the system. 

(ii) The sudden change regime without frozen discord displayed in Fig. [T] (b) . 
At first both classical and quantum correlations decrease till, at a sudden change 
point (emphasized in the figure by the first arrow on the left), classical correlation 
becomes constant while quantum correlation exhibits a discontinuous change in 
the amplitude of the damped oscillations. At the following sudden change point, 
indicated by the second arrow, classical correlation starts again to oscillate and 
quantum discord changes back to its previous rate of oscillation. 

We mention in passing that, for all the regions of parameters, entanglement 
exhibits damped oscillations, similar to those of the quantum mutual information, 
eventually showing sudden death 

4. Geometrical interpretation of frozen discord 

In Ref. Modi et al. propose a new definition for the correlations of a quantum 
system. Their idea is to put all the types of correlations on the same footing defining 
correlations in terms of distances (measured by relative entropy) between states. 
In analogy to the relative entropy of entanglement, quantum discord is defined as 
the distance of the state under study to its closest classical state, and classical 
correlation is the distance from such a classical state to its closest product state. 

In general the definition by Modi et al. does not coincide with the original 
operational one, however, for Bell-diagonal state the two formulations lead to the 
same results. Therefore we can interpret our results in the light of the correlations 
as distance measure approac h. 

As already noted in Ref. ISU, given a Bell-diagonal state p = X]i=i<i 
with the Bell states and Ai, (taken in decreasing order) the weights of each 
component, the sudden transition between classical and quantum decoherence takes 
place when the weight A3 of component $3 becomes larger than the weight A2 of 
component <I>2. This can be understood noting that the expression for the closest 
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classical state to Bell-diagonal states is 



1 - q{t) 



(19) 



2 



1=1,2 



where q{t) — Ai(t)+A2(t). Therefore when A3 becomes bigger than A2 the expression 
of the closest classical state changes abruptly. It is actually of interest to ask where 
those classical states are located in the Hilbert space of the system, and whether 
such a class of classical states display a particular geometry. 

To fix the ideas we first focus on the Markovian case illustrated in Fig. [2] (a). 
There the solid black line represents the trajectory of the state p{t) in a schematized 
Hilbert space, and the dotted red line represents the trajectory drawn by the closest 
classical state X^'^(t) (where CD stands for constant discord). In terms of the 
distance (relative entropy or trace distance), this trajectory is parallel to the one 
traced by p{t) meaning that the discord between the two lines is constant. The green 
square on the p{t) trajectory indicates the state psT at which the sudden transition 
takes place, in correspondence of A2 = A3. The state psT has two different classical 
states with equal distance: one is the last state at the end of the red dotted line, 
i. e., X^'°(Oi the other one is indicated by the red sphere at the right end of the 
black line X^^(0- Since p{t) keeps traveling the black line from left to right, after 
the transition (t) (DD meaning decaying discord) becomes the closest classical 
state. Interestingly x^^(^) is also the asymptotic state of the dynamics and does 
not evolve in time. 

The non-Markovian case is displayed in Fig. [2] (b) , there the meaning of lines 
and symbols, and the structure of the set of closest classical states remain the same 
as in Fig. [2] (a) but the time evolution of p{t) changes. Due to the non-Markovian 
memory effects, the state of the system oscillates around the transition point and 
passes many times from both of the directions the transition state psT- After the 
first passage of p{t) through psr from left to right, p{t) hits in a finite time X^^(i)) 
which in the Markovian case was reached only in the asymptotic limit. Now the 
memory of the reservoir comes into play and the state travels back part of the 
trajectory previously drawn, this is pictured by the thick dashed blue line partially 
overlapping with the black one. This time the state psr is passed from right to left, 
so p{t) enters again in the constant discord region. Eventually the direction of the 
trajectory of p{t) in the Hilbert space changes again and the state enters again the 
changing-discord-region and remaining there. 

5. Conclusive Remarks 

In this paper we have studied the dynamics of quantum and classical correlations 
for a system of two qubits locally interacting with non-Markovian non-dissipative 
environments. To solve the dynamics of the single qubits we extended the model by 
DafFer et at to the two-qubit case, when each qubit is subjected to a local colored 
noise dephasing channel. For initial Bell-diagonal states the dynamics of correla- 
tions exhibits, three different dynamical regimes: (i) constant classical correlations; 



June 22, 2011 0:32 WSPC/INSTRUCTION FILE nonMarklJQI 



Frozen discord in non- Markovian dephasing channels 11 

(ii) damped oscillations of quantum and classical correlations (iii) sudden change 
behavior with possibility for frozen discord including multiple sudden transitions 
of decoherence. The memory effects of the non-Markovian environment bring into 
play non-trivial dynamical features and allow multiple transition points between 
the classical and quantum decoherence. 

We have also proposed a geometrical interpretation of the phenomenon of frozen 
discord in terms of trajectory of the state under investigation and its closest classical 
state. In the transition point between the classical and quantum decoherence, the 
location of the closest classical state changes in discontinuous way. In the Markovian 
case, this transition point is crossed only once, where as in the non-Markovian case, 
the memory effects allow several crossings to both directions. 

Very recently Jin-Shi et al.l^have observed the dynamics of the correlations of 
two photons in a non-Markovian dephasing environment and their results display 
the same features as described here. In general, we think that result s fo r dephas- 
ing reservoirs shown here and the corresponding experimental results present a 
considerable step towards more comprehensive understanding of the dynamics of 
quantum correlations in the open system context. 
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Fig. 2. (Color online) Schematization of the trajectories of both the state of the system and its 
closest classical state in the Hilbert space in the (a) Markovian, (b) non-Markovian case. 
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